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! Abstract 



A quantum perfect lattice action in four dimensions can be derived analyt- 



' ically as a renormalized trajectory when we perform a block spin transfor- 

ms ' 

mation of monopole currents in a simple but non-trivial case of quadratic 
' monopole interactions. The spectrum of the lattice theory is identical to that 

o' 

<^ ' of the continuum theory. The perfect monopole action is transformed exactly 



into a lattice action of a string model. A perfect operator evaluating a static 



potential between electric charges is also derived explicitly. If the monopole 



interactions are weak as in the case of infrared SU (2) QCD, the string interac- 



^ ' tions become strong. The static potential and the string tension is estimated 



analytically by the use of the strong coupling expansion and the continuum 
rotational invariance is restored completely. 
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To obtain the continuum limit is crucial in the framework of lattice field theories. A 
block spin transformation which is one of renormalization group transformations is adopted 
usually on the lattice. A quantum perfect lattice action is an action on the renormalized 
trajectory on which one can take the continuum limit. 

In principle, we obtain the renormalized trajectory when we perform infinite steps (the 
n — > oo limit) of block spin transformations for fixed physical length b = na where a is the 
lattice constant. But this is actually impossible in ordinary cases. What we can do in actual 
simulations is to approach the renormalized trajectory carrying out as many steps of block 
spin transformations as possible on a finite lattice A^ 4 . If the effective action S(n, a, N) 
obtained satisfies well the two conditions that (1) S(n,a,N) is a function of b = na alone 
(the scaling behavior and volume independence) and (2) the continuum rotational invariance 
is satisfied, then the effective action could be regarded as a good approximation of the 
renormalized trajectory. The first condition can be checked when we compare S(n, a, N) 
themselves for various a, n and A^ 4 . But to test the rotational invariance, we have to 
determine the correct form of physical operators on the blocked lattice. It is the perfect 
lattice operator on the renormalized trajectory which reproduces the continuum rotational 
invariance. To find the perfect lattice operator is highly nontrivial, too. 

The purpose of this note is to give a simple but a non-trivial lattice model composed of 
general monopole quadratic interactions alone with which a block spin transformation can 
be done analytically. The renormalized trajectory and the perfect operator corresponding 
to a potential between static electric charges can be derived analytically. This is similar 
to the blocking from the continuum theory as developed by Bietenholz and Wiese Jl[. The 
spectrum of the lattice theory is the same as in the continuum theory. The continuum 
rotational invariance is shown exactly with the operator. In addition, this model is very 
interesting, since the effective monopole action obtained after an abelian projection of pure 
SU(2) lattice QCD is known to be well dominated by such quadratic monopole interactions 
alone in the infrared region [0-0] . 

Let us start from the following action composed of quadratic interactions between mag- 
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netic monopole currents. It is formulated on an infinite lattice with very small lattice 
constant a: 



S[k] = E k^Dois-s')^'). 



(1) 



Since we are starting from the region very near to the continuum limit, it is natural to assume 
the direction independence of D (s — s'). Also we have adopted only parallel interactions, 
since we can avoid perpendicular interactions from short-distance terms using the current 
conservation. Moreover, for simplicity, we adopt only the first three Laurent expansions, 
i.e., Coulomb, self and nearest-neighbor interactions. Explicitly, Dq(s — s') is expressed as 
/3A^ 1 (s — s') + a5 SyS < + 7 A L (s — s') . Here A L (s — s') = — ^ d^d' 5 S>S / and d(d') is the forward 
(backward) difference. Including more complicated quadratic interactions is not difficult. 

How to evaluate a static potential between electrically charged particles is a problem. 
It is known that the theory with the above action (TJ) is equivalent to an abelian gauge 
theory of the Villain form |||7|]. In this model, it is natural to use an abelian Wilson loop 
W(C) = exp i J2c(@n(, s )i -^( s ))> wnere $Ai( s ) is an abelian angle variable of the modified 
Villain action. Also the theory with the above action (fl|) can be rewritten |7[ in the lattice 
form of the modified London limit of the dual abelian Higgs model [Q. The static potential 
is evaluated by a 'tHooft operator in the model. However the expectation values of both 
operators are not completely equivalent, although the term of the area law is the same 0. 
When use is made of BKT transformation |T0|jrT| ,[7|j9[1 , we see that the area law term is given 
correctly also by the following operator in the monopole action []: 



Using the monopole definition a la DeGrand-Toussaint, we can prove it also directly from abelian 




(2) 



N,(s, Sj) = £ A L 1 ( S - s') - e^daS^ + A), 



(3) 



s' 



Wilson loops IHl. 
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where S^(s' + p) is a plaquette variable satisfying d'^SjJs) = J 7 (s) and the coordinate 
displacement fi is due to the interaction between dual variables. It is possible to prove that 
any choice of S^(s) for fixed electric currents J-y(s) gives the same value in the continuum 
limit a — > 0, since the difference is a closed surface and the exponent of W m (C) for the 
difference is just the four- dimensional linking number times 2iri. Hence we take a flat 
surface for Si (s) in the following. Since the area law term is the same, let us consider only 
the operator (fj) in the following. The details of the definition of the operator evaluating 
the static potential are discussed in Ref. ||. 



Now let us define a blocked monopole current [14 



n-l 



i,j,l=0 

= %(s (n) )- (4) 

With this definition, the current K^(s^) on the coarse lattice with a lattice distance b = na 
satisfies the current conservation d'^K^(s^) = Y^^K^i 3 ^) ~ K^s^ — bfi)) = 0. The block 
spin transformation is expressed as 



Z[K, J}= £ exp - K( s ) D o(s ~ s')M s ') + 2 ™ E ^(s)k^s) 

fc M = -oo I, S,s',fl S,fl ) 

e;fc M =o 

x5(K,(s^)-B kfl (s^)y (5) 

The vacuum expectation value of the Wilson loop (0) is written in terms of K^(s^) as 
follows: 

oo I oo 

(W m (C))= E z l K > J \ / E^°l- ( 6 ) 

8^=0 S;A' M =0 

Introducing auxiliary field and 7, we rewrite the constraints d'^k^ = and K fM (s^) = 
Bk As^ n '). Then we change the integral region of 7 and <p from the first Brillouin zone to the 
infinite region, since the monopole currents take integer values. Making use of the Poisson 
sum rule and recovering dimensional lattice constants a and h = na, we get 
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Z[K,J] oc r V 1 exp\tb i J2 l^bs^K^bs^)) 

00 roo roo 1' 

x Z / V $ \ VF ex n -a 8 J2 F ^ as ) D o( as - as^F^as') 

aiez 

+ia 4 J2[^N^(as) + <P(as)d'^ + 2ivl^as)] F^as) - a 4 £ n^(nas^)B Ffl (s^)\ . 



S,H s (n) 



(7) 



Since we take the continuum limit a — > finally, I = alone may remains in the sum with 
respect to al e Z. Carrying out explicitly integrals with respect to F, <p, 7 and taking the 
continuum limit a — > 0, we obtain the expectation value of the operator and the effective 
action on the coarse lattice: 



(W(C)) = expl-vr 2 rd A xd 4 yY,N,(x)D Q \x - y)N lx (y) 

+tt 2 6 8 Bp(bsW)A% F -\bsW -bsW)B v (bsW) 



s (n), s (n)' 



x 2 expi -6 8 J2 K^bs^A^J' 1 ^ - bs {n) ')K u (bs {n) ') 

6»if M (6s) = -oo I s (n) 3 (n)' 

+2irib 8 £ B^bsW)A% F -\bsW-bsW)K v (bsW)\ I £ Z[tf,0], (8) 

s (n) ]S (n)' J/ i,3_Ff M (bs) = -oo 



where 



B„(bs {n) ) = lim a 8 J] n^(6s (n) - as)A Mi ,(as - as')N v (as'), (9) 

n^oo S,s',U 

n^(6s n - as) = —5 (nastf* + (n - I) a - as M ) x [J ("asl* + /a - as t ) , (10) 

n ' i&n) \i=o J 

A^as - as') = |^ - ^^7} D \as - as'). (11) 
A^ F_1 (&s(") — bs^') is a gauge-fixed inverse of the following operator: 

A'^(bs (n) - fe (n) ') = a 8 ]T n^(6s (n) - as)A M ,(as - as')lM&s (n) - as'). (12) 

s,s' 

Here we have used Y.^d^N^(s) = and have adopted a gauge including \{d^ / y IJ ,(bs^)} 2 in 
the integral with respect to 7. 
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The momentum representation of the gauge fixed propagator is given explicitly by 



(13) 



where p^ = 2sin(p^/2) and A' (p) is written as 



4*(P)=(II £ ){D 1 (p + 2nl) 

\i=l l i= —oo) \ 

(ni^) 2 



(p + 27r/) M (p + 27r/) t 
E,(p + 2vr/)2 



(p + 27rl)^p + 2nl) u 
ni(p + 27r0? 



(14) 



From the gauge invariance condition J2fid'^A'^ v (s — s') = 0, we get ^Z^p^A' (p) = 0. The 
inverse of A'^J{p) is as follows: 



AfJ-\p) = D, v {p) + ^^-^e 



1 PyPv i( Pu - Pv )/2 

A (EiPi) 2 



(15) 



where D^p) is the A independent part of A^, F_1 (p)). To derive the explicit form of D^ v (p) 
is not so easy. However, evaluating the determinant and cofactors of the matrix A'^(p), we 
can get 

3 



D 



hi/ 



x 



A' aa ) 3 — 3 (J2 a A' aa ) (J2ab A' ab A' ba ) + 2 (X) a &c ^a&^ic^4a)} 

^la ) ~~ y?j A' ab A' ba J > 5^ — 2 ( A' aa j + 2 A'^ ia A m 

. a J \ a b /I \ a / a 



- £4- 



(16) 



Since ^^d 1 K^ibs 1 - 71 ') = ^2^8' B^bs^) = 0, Eq.(|8|) is independent of the gauge parameter 
A. 

Now let us evaluate the spectrum of the monopole current K^(s^) on the coarse lattice. 
Define an operator with definite spatial momentum p: 

K i(Plx 4 = J ^Ki(p,p 4 )e 
Then the correlation function is written as 

(^(x,0)^(£L 4 > 



-n (2tt) 4 7- w 2tt 



(17) 
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Since the monopole action on the b-lattice is written as K lx {s)D lxv {s — s')K u (s'), we see 
that the spectrum is essentially fixed by the gauge invariant part of the inverse of D^^s — s'). 
In (|T7D, it is 5 A {k + p)A' ii (p, p 4 ), where 



II E ^(p + SttZ) 
(p + 27r0i(p + 27r0<l 



E,(p+2tt/) 



IL(p + 27rZ) 



The integral 



2tt 



can be performed when we change p 4 + 27r/ 4 — > p 4 and J2u I-t+2nh ~ * /-^oo- Here 



D -» = ft 



mt 



rrin 



p 2 + m 2 p 2 + m 2 y ' 



where k, m± and m-2 are expressed by the original couplings in Eq.(|T]) as K[m\ — m 2 ) = 7 1 , 
m\ + m\ — a/7 and m\m\ = f3/ r y. Hence the p 4 integral gives us 



where + 27r/) 2 = — p\ — (p + 2irl) 2 + mf. The spectrum is identical with the one of the 
continuum. 

The above monopole action can be transformed exactly into that of the string model 
T3 , l6] . When use is made of the Poisson sum rule, we write the monopole part of Eq.(|]) as 



s,s' 



+00 



K fJi (s)=—oo \ s,s' 

+i'£F fM (s)\d^(s) + 1*Y, D A S - s')B v (s') 



S.fl 



S' ,fl,U 



f +TT 00 fir 

/ V<f>{s) E ex P - 7 E [ d M s ) + 2vr/ M (s) + 2nJ2 D^{s - s 1 )B a (s l ) 

w « M (s)=-oo I 4 s,s> si,n,a 

D~}{s - s') [d v( j}{s') + 2irl„{s') +2nJ2 - ^)Bp{s 2 ) 



(18) 



Performing BKT transformation and Hodge decomposition, we obtain 

ln(s) = s M (s) + <V( S ) 

= d,{ ~ E ^l\s'd'Ms') + r,(s>) } + E ^AI^OO, (19) 

s' a' 

where cr W(U (s) = c^s^i is the closed string variable satisfying the conservation rule 

<9[a0>] = <9a<V + ^ct^q + d v a ail = 0. (20) 

The compact field 0(s) is absorbed into a non-compact field <Pnc{ s )- Integrating out the 
auxiliary non-compact field, we see 

oo ( 

(18) = £ exp - tt 2 £ ct^c^LT 1 ^ - Sl )A L 2 ( Sl - s'K^s') 



°"/Jl/(s) = — CXD 
S [ C « CT , 1 >'] (S) = 



-2tt 2 E^(s)^ A Z 1 ( s - s')^(s') - ^E^OO^W* - s')^(s')}- (21) 

The term independent of the string variable exactly cancels the second classical term of 
Eq.©. We find finally 

(W m (C)) = ^expl-n 2 J°£xd i yJ2N fl (x)D 1 (x-y)N ll (y)\ 



oo 

X 



E exp <^ - tt 2 E ^Md^D^s - s i) A L 2 (*i " s'K/^O 



0/j,i/(s) = — OO V S f s' 



s.s' ) 



(22) 



It is very interesting that we can evaluate analytically the potential between the static 
electric charges when the monopole action on the dual lattice is in the weak coupling region 
for large b as realized in the infrared region of pure SU(2) and SU(3) QCD. Then the string 
model on the original lattice is in the strong coupling region. As shown later explicitly, the 
potential between the static electric charges is then evaluated mainly by the first classical 
part of Eq.(^) alone. Hence let us evaluate first the classical part. Since the classical part 
is written in the continuum form, the continuum rotational invariance for any b lattice site 
is trivial. 
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The plaquette variable S a p in Eq.(|D for the static potential V(Ib, 0, 0) is expressed by 



Sap{z) = 8 al 8 p4 8(z 2 )8{z 3 )9(z 1 )9(Ib - Zl )9(z 4 )9{Tb - z 4 ). 

Also the variable S a p for the static potential V(Ib, lb, 0) is given by 

S a p{z) = (6 al 6p4 + 8 a2 8 p ^j8(z ? )6(z i )6(Tb - z 4 ) 

x9(zx)9(Ib - z 1 )9(z 2 )9(Ib - z 2 )8{z 1 - z 2 ). 



(23) 



(24) 



Let us evaluate V(Ib, 0, 0) as an example. The Fourier transform of S a /3(z) in this case 



is 



rib r Tb 

S aP {p) = (8 aA 8 pi -8 al 8^) I dz 1 e-* n " 1 / dz 4 e^ z \ 



/ \ ( 2 \ - pi bl 

[8046/31 - 8 al 8f3 4 ) I — I e 1 2 sin 



vibi . ,p\Ib. ( 2 \ .-MM" . ,PiTb. 



(25) 



Since we study large T and large b behaviors, we use the following formula: 

lim (— — — — ^ = 7iT8(a). 

t^oc \ a J 



(26) 



We get 



(W{Ib,0,0,Tb)) — > exp|-7r 2 (T/6 2 ) f 

T 1 — k^v^ J 



T^oo 
b— *oo 



d 2 p 



ADr 



(0,P2,P3,0) 



(27) 



Similarly we can evaluate (W(Ib, Ib,0,Tb)) from the classical term. The static potentials 
V(Ib, 0, 0) and V(Ib, lb, 0) can be written as 

d 2 p r 1 



(2tt) : 



7(76,0,0) =7i 2 {Ib) J 



nnlb mi 
= In — , 

2 m 2 

7/6,76,0 =^ In—. 

2 m 2 



(0,P2,P3,0), 



(28) 
(29) 



The potentials from the classical part take only the linear form and the rotational invariance 
is recovered completely even for the nearest 1 = 1 sites. The string tension from the classical 
part is evaluated as 
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This is consistent with the analytical results || in Type-2 superconductor. The two constants 
mi and m2 may be regarded as the coherence and the penetration lengths. 

Next let us evaluate the quantum fluctuation coming from the interaction of the string 
variable and the classical source. Since we have introduced the source term corresponding 
to the Wilson loop on the fine a lattice, the recovery of the rotational invariance of the static 
potential is naturally expected also for the quantum fluctuation. Hence here we evaluate the 
quantum fluctuation for the flat Wilson loop W(Ib, 0, 0, Tb). Then it is to be emphasized 
that the same static potential for the flat Wilson loop can be obtained for /, T — > oo when 
we consider the naive Wilson loop operator on the course b lattice instead of that on the 
fine lattice @: 

W m (C) = exp Olm N^s, S J )K^s)) , (31) 
N^Sj) = ^A L \s - s') l -e^d a S J M {s' + fi), (32) 

s 1 " 

where SiAs' + fi) is a flat plaquette variable satisfying dpS]L(s) = J y (s) and J 7 (s) is the 
electric current on the course lattice. 

Similar arguments as above shows that the expectation value of the operator W m (C) is 
expressed as follows: 

(W m {C)) = | exp | - ^J2Y,N»(s)D; u \s - s')N v (s>)} 

OO ( 

x ex P 1 - ^ 12 a ^(s)d a d' p D-^(s - S!)A^ 2 (si - s')a u/3 (s') 



CTflu(s) — — 00 



-27r 2 ^a^( S )^A^( S - s x )D va { Sl - s')N a {s')\. (33) 

s,s> J 

Here we note that D~^(s — s') is given by Eq. (|I2"D and that the Fourier transform Sp^k) of 
S"^ 7 (s) is similar to Eq.(f25|) with the momentum defined on the course lattice. Hence the 



first classical part becomes 
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TV 2 



J2N,(s)D^(s-s>)N„(s') 



x S , i 4 (A;)S , i4(-A;) 2_, + 2tt/) 2 JA>(fc + 2 ^) 

x IIi# M 7r < (A; + 27rO<(ife + 27rO, (34) 



where 



*<*) = 5^-**. (35) 

Changing the integral variable as + 27r/ — > k, we can absorb the summation with respect 
to I using the integral over the infinite momentum range. When we use Eq.(p6|) for large / 
and T, we find the classical part Q53j ) agrees exactly with (|2T|). 
Similarly 

a> 



-7T (2tt) 



/ (k + 2nl),(k + 2nl)A 
£(V )D (k + 2nl) 



x U^7r*(k + 2irl)irj(k + 2Til)e l ^2ie yal ^l\k) sin e ^ , > e » ( 36 ) 

y 2 J TT\\k)Tii\k) 



where 



„ m = sm^J/2) e _, fciJ/2) ^ (Jfe) = sin(fc 4 T/2) e _ ifc4T/2 _ (3?) 



When use is made of 

sm sin — = sin — sin , 38 

V 2 J \2j { 4 J { 4 ) 1 J 

we get for large / and T 

7rJ(ife)7fi(A;)7r:(A;)7f4(Jfe) = (27r) 2 5(A; 1 )5(fc4)e i(1 - /)fcl/2 e i(1 - T)fc4/2 . (39) 
We find for large I and T that £) a / D~^(s — s')N u (s') is equivalent to B^s) in (g). 
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Hence the expectation value of the naive Wilson loop (W m (C)) in (|33"D coincides with 



that of the perfect operator in (p2|). Now we introduce the dual string variable V as follows: 



<V(s) = ^va^CTapis + & + $), 

*a a/3 (s) = *a a p + S a p(s). (40) 



Then the expectation value can be expressed simply as 



where 



5(*<t) = vr 2 ^(fiaAZ 1 ^ - s^-e^a^s, + a + ft))D^(s - s') (42) 



s,s' 



x^Aj/V - 82)^^5* a v5 {s 2 + $ + v)). (43) 

The strong coupling expansion can be shown in Fig.|I|. The leading term is the same as the 
classical contribution in (0). The next-leading term is a house-type diagram with one lxl 
cube attached on the flat surface. Then the open surface variable *a has four more plaquettes 
than the leading one. If the self coupling term between the string variables *a is dominant 
as in SU{2) QCD, the next to leading term is estimated as exp(—<j d ITb 2 — 4I1(0)6 2 ), where 
a c i is the string tension from the classical part fl3"0|) and 11(0) is the self coupling constant. 
Considering the entropy factor AIT, we get 

(W m (C)} = exp{-a cl ITb 2 } + AIT exp{-a cl ITb 2 - AU(0)b 2 } + ■ - 
~ exp{-a c/ /T6 2 } exp{4/T e - 4n(0)fc2 }, 

= exp{-aITb 2 }. (44) 



Hence the string tension becomes 



a = a d -le-W 2 . (45) 



Applications to actual pure 577(2) and 577(3) QCD will be published elsewhere PJT8 
and the quantum fluctuation term will be found to be very small there. 
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FIGURES 



FIG. 1. The strong coupling expansion of the Wilson loop calculation. 
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